APS/V2-15/10/02 



o 
o 

(N 
< 



0) 

o 
Or 

r/3 ' 
1/3 ■ 

o ■ 

C/3 . 
>v 
43 ■ 
Oh! 



> 
o 

OO 

o 



(N 

o 

t/3 

O 

C/3 1 

>>: 
43 : 

Oh' 



X 



Intermittency, scaling and the Fokker-Planck approach to 
fluctuations of the solar wind bulk plasma parameters as seen by 

the WIND spacecraft. 

Bogdan HnatQ Sandra C. Chapman, and George Rowlands 
Physics Department, University of Warwick, Coventry, CV4 7AL, UK. 
(Dated: February 2, 2008) 

The solar wind provides a natural laboratory for observations of MHD turbulence over extended 
temporal scales. Here, we apply a model independent method of differencing and rescaling to 
identify self-similarity in the Probability Density Functions (PDF) of fluctuations in solar wind bulk 
plasma parameters as seen by the WIND spacecraft. Whereas the fluctuations of speed v and IMF 
magnitude B are multi-fractal, we find that the fluctuations in the ion density p, energy densities 
B 2 and pv 2 as well as MHD-approximated Poynting flux vB 2 are mono-scaling on the timescales up 
to 26 hours. The single curve, which we find to describe the fluctuations PDF of all these quantities 
up to this timescale, is non-Gaussian. We model this PDF with two approaches- Fokker-Planck, 
for which we derive the transport coefficients and associated Langevin equation, and the Castaing 
distribution that arises from a model for the intermittent turbulent cascade. 



PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Statistical properties of velocity field fluctuations 
recorded in wind tunnels and these obtained from solar 
wind observations exhibit striking similarities 
unifying feature found in these fluctuations is fractal or 
multi-fractal scaling. The Probability Density Function 
(PDF), unlike power spectra that do not reveal intermit- 
tency, show a clear departure from the Normal distribu- 
tion when we consider the difference in velocity on small 
spatial scales yl 13 while large scale features appear to 
be uncorrelated and converge toward a Gaussian distri- 
bution. These similarities suggest a common origin of 
the fluctuations in a turbulent fluid and the solar wind. 
The approach is then to treat the solar wind as an active 
highly nonlinear system with fluctuations arising in situ 
in a manner similar to that of hydrodynamic turbulence 

Kolmogorov's K41 turbulence theory was based on the 
hypothesis that the energy is transferred in the spec- 
tral domain at a constant rate through local interaction 
within the inertial range. This energy cascade is self- 
similar due to the lack of any characteristic spatial scale 
within the inertial range itself. These assumptions led 
Kolmogorov to his scaling law for the moments of veloc- 
ity structure functions 5 1 " = (\v(r + £) — v(r)\ n ) oc 
(e£)™/ 3 , where n is the ra-th moment, £ is a spatial scale 
and e represents energy transfer rate. Experimental re- 
sults do not confirm this scaling, however, and modifi- 
cations to the theory include intermittency '9j by means 
of a randomly varying energy transfer rate e. In this 
context, empirical models have been widely used to ap- 
proximate the shapes of fluctuation PDFs of data from 
wind tunnels ^(j as we U as the solar wind; see for ex- 
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ample 0, ■ The picture of turbulence emerging from 
these models is much more complex then has been sug- 
gested by the original Kolmogorov theory. It requires a 
multi-fractal phenomenology to be invoked as the self- 
similarity of the cascade is broken by the introduction of 
the intermittency. 

Recently, however, a new approach has emerged where 
the presence of intermittency in the system coincides with 
statistical self-similarity, rather than multi-fractality, in 
the fluctuations of selected quantities; these also exhibit 
leptokurtic PDFs. An ex amp le of this statistical inter- 
mittency was discussed in [l8l , where a Levy distribution 
was successfully fitted to the fluctuation PDFs of the 
price index over the entire range of data. Such a distribu- 
tion arises from the statistically self-similar Levy process 
also characterized by enhanced (when compared with 
a Gaussian) probability of large events. Recently Q 
reported similar self-similarity derived from the scaling 
of the solar wind Interplanetary Magnetic Field (IMF) 
energy density fluctuations calculated from the WIND 
spacecraft dataset. Here, we apply a model-independent 
and generic PDF rescaling technique to extract the scal- 
ing properties of the solar wind fluctuations directly from 
the data. The aim is to determine a set of plasma param- 
eters that exhibit statistical self-similarity and to verify 
the nature of the PDF for their fluctuations. We con- 
sider the following bulk plasma parameters: magnetic 
field magnitude B, velocity magnitude v, ion density p, 
kinetic and magnetic energy density (pv 2 and B 2 ) and 
Poynting flux approximated by vB 2 . Such an approxi- 
mation of the Poynting flux assumes ideal MHD where 
E = v x B. We find that the PDFs of fluctuations in p, 
B 2 , pv 2 and vB 2 exhibit mono-scaling for up to 10 stan- 
dard deviations, while B and v are clearly multi-fractal 
as found previously 0, 0] . 

The mono-scaling allows us to derive a Fokker-Planck 
equation that governs the dynamics of the fluctuations' 
PDFs. The Fokker-Planck approach provides a point of 
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FIG. 1: Unsealed PDFs of the ion density fluctuations. Time 
lag r = 2 fe x 46s, where k = 0, 1, 2, .., 14. The standard devia- 
tion of the PDF increases with r. The error bars on each bin 
within the PDF are estimated assuming Gaussian statistics 
for the data within each bin. 
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FIG. 2: Scaling of the peaks P(0, r) of the PDFs for all quan- 
tities under investigation: o corresponds to SB 2 , □ ion den- 
sity dp, o kinetic energy density S(pv 2 ) and A Poynting flux 
component 8(vB 2 ). The plots have been offset vertically for 
clarity. Errors are estimated as in Fig. 



contact between the statistical approach and the dynam- 
ical features of the system. This allows us to identify the 
functional form of the space dependent diffusion coeffi- 
cient that describes the fluctuations of these quantities as 
well as to develop a diffusion model for the shape of their 
PDFs. We also consider a Castaing model where fluctua- 
tions are assumed to arise from a varying energy transfer 
rate e in the nonlinear energy cascade, with Gaussian 
distribution for ln(e). The paper is structured as fol- 
lows: in section II we will describe the dataset used for 



this study as well as the rescaling procedure. In section 
III the results of the rescaling will be presented. Two 
possible models of the fluctuations will be discussed in 
Section IV. Finally in Section V we will summarize all 
results discussed throughout this paper. 



II. DATA AND METHODS 

A. The Dataset 

The solar wind is a supersonic, super- Alfvenic flow of 
incompressible and inhomogeneous plasma. The WIND 
spacecraft orbits the Earth-Sun LI point providing a 
set of in situ plasma parameters including mag netic 
field measurements from the MFI experiment [Ig and 
the plasma parameters from the SWE instrument |17| . 
The WIND solar wind magnetic field and key parameter 
database used here comprise over 1.5 million, 46 second 
averaged samples from January 1995 to December 1998 
inclusive. The selection criteria for solar wind data is 
given by the component of the spacecraft position vector 
along the Earth-Sun line, X > 0, and the vector magni- 
tude, R > 30 RE. The data set includes intervals of both 
slow and fast speed streams. Similar to other satellite 
measurements, short gaps in the WIND data file were 
present. To minimize the errors caused by such incom- 
plete measurements we omitted any intervals where the 
gap was larger than 2%. The original data were not av- 
eraged nor detrended. The data are not sampled evenly 
but there are two dominant sampling frequencies: 1/46 
Hz and 1/92 Hz. We use sampling frequency f s of 1/46 
as our base and treat other temporal resolutions as gaps 
when the accuracy requires it (r < 92 seconds). 



B. Differencing and Rescaling Technique 

Let x(t) represent the time series of the studied sig- 
nal, in our case magnetic field magnitude B, velocity 
magnitude v, ion density p, kinetic energy density pv 2 , 
magnetic field energy density B 2 or the Poynting flux 
component approximated by vB 2 . A set of time series 
Sx(t, t) = x(t + t) — x(t) is obtained for each value of the 
non-overlapping time lag r. The PDF P(Sx, r) is then 
generated for each time series 6x(t, r). Fig. ^shows the 
set of such raw PDFs of the density fluctuations for time 
lags between 46 seconds and ~ 9 days. A generic one pa- 
rameter rescaling method 01 1S applied to these PDFs. 
We extract the scaling index a, with respect to r, directly 
from the time series of the quantity Sx. Practically, ob- 
taining the scaling exponent relies on the detection of a 
power law, P(0, r) cx r~ Q , for values of the raw PDF 
peaks and time lag r. Fig. |21 shows the peaks P(0, r) of 
the unsealed PDFs plotted versus r on log-log axes for 
the four bulk plasma parameters. We see that the peaks 
of these PDFs are well described by a power law r~ a for 
a range of t up to ~ 26 hours. We now take a to be the 
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FIG. 3: One parameter rescaling of the PDF for the fluctu- 
ations in the magnetic field energy density B 2 . The curves 
shown correspond to r between 2 minutes and 26 hours. 
Error bars as in Fig. 

scaling index and attempt to collapse all unsealed PDFs 
P(Sx,t) onto a single curve P s (Sx s ) using the following 
change of variables: 

P(Sx,t) =r- a P s {SxT- a ). (1) 

A self-similar Brownian walk with Gaussian PDFs on all 
temporal scales and index a — 1/2 is a good example 
of the process where such collapse can be observed (see 
e.g. [l8j)- For experimental data, an approximate col- 
lapse of PDFs is an indicator of a dominant self-similar 
trend in the time series, i.e., this method may not be 
sensitive enough to detect multi-fractality that could be 
present only during short time intervals. One can treat 
the identification of the scaling exponent a and, as we 
will see, the non-Gaussian nature of the rescaled PDFs 
(P s ) as a method for quantifying the intermittent char- 
acter of the time series. Another possible interpretation 
of the rescaling is to treat P(Sx, r) as the self-similar 
solution of the equation describing the PDF dynamics. 
The mono-scaling of the fluctuations PDF, together with 
the finite value of the samples' variance, indicates that a 
Fokker-Planck approach can be used to express the dy- 
namics of the unsealed PDF in time and with respect 
to the coordinate Sx ^J- I n section IV we will use the 
Fokker-Planck equation to develop a dynamical model 
for the fluctuations observed in the solar wind. 

Ideally, we use the peaks of the PDFs to obtain the 
scaling exponent a, as the peaks are statistically the most 
accurate parts of the distributions. In certain cases, how- 
ever, the peaks may not be the optimal statistical mea- 
sure for obtaining the scaling index. For example, the 
B z component of the solar wind magnetic field is mea- 
sured with an absolute accuracy of typically about 0.1 
nT. Such discreteness in the time series introduces large 
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FIG. 4: As in Fig.[3]for ion density fluctuations 8 p. 
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FIG. 5: As in Fig. for kinetic energy density fluctuations 
5(pv 2 ). 

errors in the estimation of the peak values P(0, r) and 
may not give a correct scaling. However, if the PDFs 
rescale, we can in principle obtain the scaling exponent 
from any point on the curve. We will illustrate this in 
the next section where we obtain the rescaling index a 
from two points on the curve P(0,t) and P(<j,t). 



III. PDF RESCALING RESULTS 

We are now ready to present results of the rescaling 
procedure as applied to the solar wind bulk plasma pa- 
rameters. Fig. ^ shows the unsealed (raw) PDF curves 
of the ion density data. These PDFs, like all others pre- 
sented in this section, were generated with the bin size 
decreasing linearly toward the center of the distribution 
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to improve the accuracy of the PDF for small fluctua- 
tions. Although the entire range of data was used to 
create these PDFs we truncated the plotted curves for 
\Sx\ > 10<t(t), where tr(r) is a standard deviation of the 
differenced time series for the specific time lag r. Fig. [21 
then shows P(0, r) plotted versus r on log-log axes for 
5x = S(p), 6(pv 2 ), S(B 2 ) and 8(vB 2 ). Straight lines on 
such a plot suggest that the rescaling Q holds at least 
for the peaks of the distributions. In Fig. [2J lines were 
fitted with R 2 goodness of fit for the range of r between 
2 minutes and 26 hours, omitting points corresponding 
to the first two temporal scales as in these cases the sharp 
peaks of the PDFs can not be well resolved. The lines 
suggest self-similarity persists up to intervals of r w 26 
hours. The slopes of these lines yield the exponents a 
and these are summarized in Table [I] along with the val- 
ues obtained from analogous plots of P(<j(t), t) versus r 
which show the same scale break and the same scaling 
exponent for S(p), S(pv 2 ), 8(B 2 ) and S(vB 2 ), to within 
the estimated statistical error. Within this scaling range 
we now attempt to collapse each corresponding unsealed 
PDF onto a single master curve using the scaling JQ. 
Figs. 15151 show the result of the one parameter rescaling 
applied to this unsealed PDF of fluctuations in p, pv 2 , 
B 2 and vB 2 respectively, for temporal scales up to ~ 26 
hours. We see that the rescaling procedure ^ using the 
value of the exponent a of the peaks P(0, t) shown in 
Fig. [2 gives good collapse of each curve onto a single 
common functional form for the entire range of the data. 
These rescaled PDFs are leptokurtic rather than Gaus- 
sian and are thus strongly suggestive of an underlying 
nonlinear process. 

The fluctuations PDFs for all mono-scaling quantities 
investigated here are nearly symmetric. This is in sharp 
contrast with the strong asymmetry of the PDF of ve- 
locity fluctuations in hydrodynamic turbulence reported 
previously in [ToL l2fj| . This asymmetry of the statistics 
for the velocity increments coincides with the highly in- 
termittent character of the flow and multi-fractal scaling 
of these fluctuations. We applied zeroth-order correlation 
functions, defined separately for the positive and the neg- 
ative branch of the PDF [2(J , to quantify the asymmetry 
of fluctuations PDFs for the solar wind. This analysis 
was performed using PDFs generated for r w 12 minutes 
(that is, within the scaling region). In the case of veloc- 



Quantity 


a from P(0, r) 


a from P(a,r) 


Approx 

Tmax 


PDF 

scales 


SB 


-0.47 ±0.02 


-0.23 ± 0.05 


26 hrs 


No 


Sv 


-0.52 ±0.05 


-0.21 ±0.06 


26 hrs 


No 


sm 


-0.43 ± 0.03 


-0.39 ± 0.08 


26 hrs 


Yes 


S(P) 


-0.39 ±0.03 


-0.37 ±0.05 


26 hrs 


Yes 


S(pv') 


-0.41 ±0.03 


-0.35 ± 0.05 


26 hrs 


Yes 


S(vB' J ) 


-0.42 ± 0.02 


-0.39 ± 0.06 


26 hrs 


Yes 



TABLE I: Scaling indices derived from P(0, r) and P(<t, t) 
power laws. 
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FIG. 6: As in Fig.Elfor Poynting flux 8{vB 2 ). 



ity increments we find that the negative moment is, on 
average, 11% lower compared to the positive one. On the 
other hand, the quantities that we have found with self- 
similar increments do not have appreciable asymmetry. 
These show differences between negative and positive mo- 
ments of about 2 — 3%, which is however, well above the 
statistical error of this procedure. 

It has been reported previously ^(j that the PDFs 
obtained from hydrodynamic turbulence have exponen- 
tial tails. These would look linear on the linear-log plots 
that are used in this paper. In the case of solar wind 
bulk plasma parameters we do not find such a clear expo- 
nential cutoff region but rather see stretched exponential 
tails of the form P(|fe|) ~ exp(— ^l&rl^). This is illus- 
trated in Fig. [7| where we plot log(log(P(<5:r))) against 
\og(5x) for all positive fluctuations of the mono-scaling 
quantities. It can be seen that, as we move away from 
the peak, these curves converge to lines and good fits can 
be obtained in the interval [2a, 10er], where a stands for 
the standard deviation. 

We can now directly compare the functional form of 
these rescaled PDFs by normalizing the curves and over- 
lying them on the single plot for a particular r within 
the scaling range. Fig. [S] shows these normalized PDFs 
P s (6x S: t) for Sx s = S(p) s , S(B 2 ) S , 5{pv 2 ) s , 5{vB 2 ) s and 
t w 1 hour overlaid on a single plot. The Sx s variable 
has been normalized to the rescaled standard deviation 
(j s (t w Ihr) of P s and the values of the PDF has been 
modified to keep probability constant in each case to fa- 
cilitate this comparison. These normalized PDFs have 
remarkably similar functional form suggesting a shared 
process responsible for fluctuations in these four plasma 
parameters on temporal scales up to T max ss 26 hours. 

It has been found previously 0] that the magnetic 
field magnitude fluctuations are not self-similar but 
rather multi-fractal. For such processes the scaling de- 
rived from P(0, r) would not be expected to rescale 
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FIG. 7: Positive tails of all self-similar PDFs for r « 30 
minutes. Solid lines show linear fits obtained in the interval 
[3cr, 6<j] and extended to [2a, 10a] fluctuations interval. 



the entire PDF. To verify this we applied the resell- 
ing procedure for magnetic field magnitude differences 
SB(t,r) = B(t + T)-B(t). Fig. EJshows the result of one 
parameter rescaling applied to the PDFs of the magnetic 
field magnitude fluctuations. We see that the scaling pro- 
cedure is satisfactory only up to ~ 3 standard deviations 
of the original sample, despite the satisfactory scaling 
obtained for the peaks P(0, r) of the PDFs (see insert 
of the Fig. |3 ■ This confirms the results of where a 
two parameter Castaing fit to values within 3 standard 
deviations of the original sample yields scaling in one pa- 
rameter and weak variation in the other. Attempts to 
improve the collapse by using information in the tails 
(values \5B\ > 3er) would introduce a significant error in 
the estimation of the scaling exponent a. We found sim- 
ilar lack of scaling in the fluctuations of the solar wind 
velocity magnitude and we show the rescalcd PDF in the 
Fig. El We stress that the log-log plots of the PDF 
peaks P(0,t) show a linear region for both velocity and 
magnetic field magnitude fluctuations (see insert in each 
figure). Their PDFs, however, do not collapse onto a sin- 
gle curve when the rescaling is applied. This lack of 
mono-scaling is evident when indices derived from P(0, r) 
and these found for P(a, r) are compared (see Tabled. 



IV. MODELLING THE DATA 

The rescaling technique applied in the previous section 
indicates that, for certain temporal scales, the PDFs of 
some bulk plasma parameters can be collapsed onto a sin- 
gle master curve. The challenge now lays in developing 
physical models that can describe the functional form of 
this curve. Here we consider two approaches. The first is 
a statistical approach where we assume that the fluctua- 
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FIG. 8: Direct comparison of the PDFs of fluctuations for all 
four quantities, o corresponds to S(B 2 ), □ ion density 8{p), o 
kinetic energy density S(pv 2 ) and A Poynting flux component 
S{vB 2 ). 



tions can be described by a stochastic Langevin equation. 
The second method is to assume the fluctuations are the 
result of the nonlinear energy cascade and derive the cor- 
responding PDF form for the rescaled PDFs (Castaing 
distribution) [Tof . 



A. Diffusion model 

The Fokker-Planck (F-P) equation provides an impor- 
tant link between statistical studies and the dynamical 
approach expressed by the Langevin equation 0. In 
the most general form F-P can be written as: 

dP 

— = V Sx (A(5x)P + B(6x)V Sx P), (2) 
or 

where P = P(8x, t) is a PDF for the differenced quantity 
8x that varies with time r, A(Sx) is the friction coeffi- 
cient and B(8x) is related to a diffusion coefficient which 
we allow to vary with dx. For certain choices of A(5x) 
and B(Sx), a class of self-similar solutions of J5J satis- 
fies the rescaling relation given by (JTJ- This scaling is 
a direct consequence of the fact that the F-P equation 
is invariant under the transformation Sx — > Sxr~ a and 
P — * Pr a . It can be shown (see Appendix A) that 
equations Q and J5J combined with power law scaling 
of the transport coefficients A(Sx) and B(Sx) lead to the 
following equation for the PDF: 



dP 
97 



d 



d{Sx) 



{5x) 1 - 1,a [a Q P + b Q 5x 



dP 
d(5x) 



, (3) 



where ao and bo are constants, a is the scaling in- 
dex derived from the data and P(5x) and 5x are un- 
sealed PDF and fluctuations respectively. Written in 
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FIG. 9: As in Fig.l^for the solar wind magnetic filed magni- 
tude fluctuations. 
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FIG. 10: As in Fig.[3]for the solar wind velocity fluctuations. 



this form equation J2J immediately allows us to identity 
the functional form of the diffusion coefficient, namely 
D(Sx) oc ((5a;) 2-1 / Q . In Appendix A we show how © 
can also be expressed as: 



bo fx . dP s 
a d(dx s ) 



P s + — (5x s )«P B = C. (4) 
a 



The partial differential equation can be solved ana- 
lytically and one arrives at the general solution in the 
form: 



P s (Sx 



f'o 



C 



-exp 



b \5x s \ a °/ b °~ r V b 
s Xb exp(^{5x' s fl a 

(Sx' s y- a °/ b ° 



(5x s fl a 
d(Sx' s )+k H(5x s ), (5) 



where ko is a constant and H(Sx s ) is the homogeneous 
solution: 



H(Sx s ) 



1 



{Sx s ) a »/ b « 



exp 



bo 



(6) 



We then attempt to fit the predicted solution to the 
normalized rescaled PDFs. The results of such a fit for 
the fluctuations of the kinetic energy density PDF is 
shown in Fig. ^] (solid line). This fit is obtained with the 
following parameters ao/bo = 2.0, bo = 10, C = 0.00152, 
ko = 0.0625 and a — 0.41 as derived from the rescaling 
procedure. We note that the figure is a semi-log plot and 
thus emphasizes the tails of the distribution - for a dif- 
ferent value of the ratio ao /bo the fit around the smallest 
fluctuations could be improved. Equation JSJ can not, 
however, properly model the smallest fluctuations as it 
diverges for Sx s — ► 0. 

Let us now assume that a Langevin equation in the 
form 



d(Sx) 
dt 



p(dx)+ 1 (dx)^t) 



(7) 



can describe the dynamics of the fluctuations. In J2J the 
random variable £(t) is assumed to be 5-correlated, i.e., 



(8) 



This condition is fulfilled in the data analysis by forming 
each time series Sx(t, r) with non-overlapping time inter- 
vals r and was also verified by computing the autocorre- 
lation function of the differenced time series. Introducing 
a new variable z = J Q X l/j(Sx')d(Sx'), equation (JJJ can 
be written as: 



dz 



7 (z) 



£(*)■ 



(9) 



One can immediately obtain a F-P e quat ion that corre- 
sponds to the Langevin equation © 19]. We can then 
compare this F-P equation with that given by 10 to ex- 
press coefficients (3{5x) and j(Sx) in terms of ao and bo 
(see Appendix B). Defining Do =< £, 2 (t) > /2 we obtain: 




and 



j(6x) 



P(Sx) = Ml-— )-ao](Sx) 1 ~' 
2a 



(10) 



(11) 



Equation Q together with definitions of its coefficients 
(|10fl and I jllj l constitutes a dynamical model for the fluc- 
tuations in the solar wind quantities. From (|10fl and 
(flip, we see that the diffusion of the PDF of fluctuations 
in the solar wind is of comparable strength to the ad- 
vection (a /&o ~ 2). We stress that the advection and 
diffusion processes that we discuss here are of the proba- 
bility in parameter space for fluctuations and do not refer 
to the integrated quantities. 
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FIG. 11: Example of the fit of the PDF functional form pre- 
dicted by a Fokker-Planck description (HJ (solid line) and a 
Castaing model (dash line) to the fluctuations PDF of the 
8{pv 2 ) bulk parameter. 

B. Castaing model 

We now, for comparison, consider a model motivated 
directly by a cascade in energy, due to Castaing. This em- 
pirical model was developed for the spatial velocity fluc- 
tuations recorded from controlled experiments in wind 
tunne ls [Tol |21l and has been applied to the solar wind 
data [llj [13 ■ The underlying idea of this approach is 
that, for constant energy transfer rate between spatial 
scales, all quantities should exhibit a Gaussian distribu- 
tion of fluctuations. The intermittency is then introduced 
to the PDF through the fluctuations of the variance a of 
that Gaussian distribution. A log-normal distribution is 
assumed for the variance a: 

QW) = ^eap d{ln{*)), (12) 

where <ro is the most probable variance of the fluctuations 
and A is the variance of Zn(cr). Combining these two hy- 
pothesis Castaing proposed the following functional form 
for the observed PDF: 

(13) 

The dashed line in the Fig. II II shows the Castaing curve 
fitted with parameters A = 1.275 and uq — 0.225 to the 
S(pv 2 ) PDF. 

We can now compare the rescaled PDFs with both 
F-P and Castaing predicted curves which are shown in 
Fig. ^2 We can see from the figure that both models 
provide an adequate fit to the S(pv 2 ) s PDF, and hence 
will also describe the PDF of other scaling bulk plasma 
parameters. Both curves, however, fall significantly be- 
low observed PDF values for |c)(/?i> 2 ) s | < 2, although the 



Castaing distribution fits the peak of the PDF reason- 
ably well (see insert in Fig. Ill|) . This departure from the 
experimental PDF, in the case of the Castaing distribu- 
tion, may reflect the difference between hydrodynamics 
and MHD turbulence. 



V. SUMMARY 

In this paper we have applied a generic PDF rescaling 
method to fluctuations in the solar wind bulk plasma pa- 
rameters. We find that, consistent with previous work, 
magnetic field and velocity magnitude fluctuations are 
multi- fractal whereas the PDFs of fluctuations in B 2 , p, 
pv 2 and vB 2 can be rescaled with just one parameter 
for temporal scales up to ~ 26 hours. The presence of 
intermittency in the plasma flow is manifested in these 
quantities simply by the leptokurtic nature of their fluc- 
tuation PDFs, which show increased probability of large 
fluctuations compared to that of the Normal distribution. 
Fluctuations on large temporal scales, r > 26 hours are 
uncorrelated in that their PDFs converge toward a Gaus- 
sian distribution. The fact that all quantities share the 
same PDF, to within errors, is also strongly suggestive 
of a single underlying process. This is also supported by 
the similar values of the scaling exponents. 

The simple scaling properties that we have found al- 
low us to develop a Fokker-Planck approach which pro- 
vides a functional form of the rescaled PDFs as well as a 
Langevin equation for the dynamics of the observed fluc- 
tuations. The model shows that both advective and dif- 
fusive terms need to be invoked to describe the dynamics 
of the fluctuations. The calculated diffusion coefficient is 
of the form D(x s ) cx {5x s ) 2 ~ 1 / a . We obtained a good fit 
of the model to our rescaled PDFs over at least 10 stan- 
dard deviations. We also examined a Castaing model for 
turbulence and found a set of fit parameters for which 
both the Castaing distribution and our diffusion model 
have nearly identical form. Since both the F-P model and 
the Castaing distribution fit our rescaled PDFs we con- 
clude that their moments should exhibit same variation 
with time lag r. 
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APPENDIX A 

Let P(8x, t) be a homogeneous function that satisfies 
scaling Our aim is to find functional form of the coef- 
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ficients A(Sx) and B(8x) for which P(Sx,t) is a solution 
of a F-P equation (J2J. Using we can now rewrite 
to read: 



where the random variable £(t) is assumed to be J- 



correlated, i.e., 



" f P.+8x. 



P s dA(Sx) A{8x) dP s 



<mat+r)>=a 2 s( T ). 



(B2) 



V d{Sx s )J t a d{Sx) t 2a d(5x s ) 

1 dB(Sx) dP s B(Sx) dP s 



t 2a d{Sx) d(Sx s ) t 3a d{5x s )' 



(Al) 



If all terms in the rhs of (|A1|) are to contribute and for 
P(Sx s ) to remain a function of Sx s only we must have: 

-i-^ = a(fa s ) and -A_Z = &(fo: s ). (A2) 

Both and B(Sx) must then be of form: 

=a (Sx) ri and = 6 (fe) I/ , (A3) 

where ao and 6o are constants. Changing variables to 
the rescaled 8x s and substituting i|A3[l into i|A2[l we ex- 
press exponents 7/ and v in terms of the rescaling index 
a derived from the data. We then obtain: 



Tj = 1 and v = 2 , 

a a 



(A4) 



which allows to write the final power law form of A{5x) 
and B(5x): 



A(5x) = ao(Sx) 



and B(6x) = bo(Sx) 



2-i 



(A5) 



Substituting these expressions into F-P equation wc 
obtain © from Section 4. Using these results the term 



d d(Sx) on * ne rns °^ f° r example, becomes: 



dA(Sx) 
d(Sx) 



1 



a (fe) 



Performing similar algebra on all terms in (jAljl we arrive 
to equation: 



d(5x s P s 



(Sx^-i ( 



d(Sx s ) d(Sx s ) 
Integrated once we obtain equation (@J 



a P s + b (Sx s 



dP s 
d(Sx s ) 
(A7) 



bo , x ■. dP s 
a d(dx s ) 



P s + —(5x s )«dP s =C, (A8) 
a 



where C is the constant of integration. 



Introducing a new variable z — J Q Sx l/j(6x')d(Sx') 1 equa- 
tion (IBlll can be written as: 



^=T(z) + Z(t), where T(z) = ®. (B3) 
dt 7(z) 



One can immediately obtain a F-P equation that corre- 
sponds to the Langevin equation l|B3|) and reads: 



dP(z,r) d , , , , . . d 2 P(z,T) , N 
■ — (T(z)P(z,r)) = D \ ( B4 ) 



dr dz 



d 2 z 



where D = er 2 /2. The probability is an invariant of the 
variable change so that P(5x)d(5x) = P(z)dz and we can 
then rewrite i|B4(l for P(Sx,t): 



d 



dP 



dr d(Sx) 



D ,(Sx) d 4?4 - P(Sx)) P + D 01 * dP 



d{8x) 



d(5x) 
(B5) 



Comparing (|B5(1 with the F-P equation © we can iden- 
tify: 



A)7 2 - {Sx^boSx, 



(B6) 



(A6) 

and then wc must demand that: 



D dj 2 (Sx) x _i 

"r^feF"^^ ^ (B7) 



In summary we have shown that the F-P equation given 
by ijsjl is equivalent to the stochastic Langevin equation 
(|7|) where coefficients (3 and 7 are given by: 



7 = 




(Sx) 1 -^ 



(B8) 



APPENDIX B 

Consider the following Langevin type of equation 
d(6x 



dt 



= /3(6x)+ 1 (Sx)at), 



(Bl) 



and 



= 



1 -2^ | - ao 



(for) 1 "* 



(B9) 
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